DESY 02-016 
SFB-CPP-04/11 
January 2004 



hep-ph/0409289 



Universal Higher Order QED Corrections 
to Polarized Lepton Scattering 



Johannes Bliimlein and Hiroyuki Kawamura 1 



Deutsches Elektronen-Synchrotron, DESY, 
Platanenallee 6, D-15738 Zeuthen, Germany 



Abstract 

We calculate the universal radiative QED corrections to polarized lepton scattering for 
general scattering cross sections in analytic form. The flavor non-singlet and singlet ra- 
diation functions are calculated to 0((aL) 5 ). The resummation of the non-singlet and 
singlet contributions to the QED-anomalous dimensions oc (aln 2 (x)) k is performed to all 
orders. Numerical results are presented for the individual radiation functions. Applications 
to polarized deeply inelastic lepton-nucleon scattering are given. 



1 Present address : KEK High Energy Accelerator Research Organization, 1-1 Oho, Tsukuba-shi, Ibaraki-ken 
305-0901, Japan. 



1 Introduction 



QED corrections to integral and differential cross sections of light charged lepton-lepton scat- 
tering or deeply inelastic lepton-nucleon scattering can be quite large in some kinematic re- 
gions [1-3]. In particular Bremsstrahlung contributions lead to significant shifts in the sub- 
system kinematics and cause large logarithmic contributions in differential distributions. The 
radiative corrections can be grouped into universal pieces, which are process-independent, and 
the remainder process dependent part. The former corrections are called universal corrections 
as they are associated with the radiating particle irrespective of the process. Two classes 
of universal corrections are known : i) the leading-logarithmic O (aL) k corrections, where 
L = ln(s/m 2 ). Here s is a typical energy scale which describes a time- or a space-like virtuality 
and m denotes the light fermion mass. This type of corrections can be obtained by solving the 
renormalization group equations for mass-factorization [4] in lowest order completely; it) the 



QED-exponentiation 



a In (x) 



k 



-terms in the all-order anomalous dimensions due to infrared 
evolution equations for initial-state radiation in lowest order [5,6] 2 . 

Unlike the case in QCD, for QED radiation the radiation sources are no extended distribu- 
tions, fi(x,Q 2 ), which fall rapidly as x — > 1. The sources are given by a 5(1 — ^-distribution 
instead. A numerical solution of the evolution equations via a Mellin transform with a complex 
argument N and a subsequent numerical inverse Mellin transform is thus not possible since the 
latter requires substantial damping as Re(iV) — > — oo, which is not the case in general here. 
To determine the solution of the respective evolution equations in QED one rather evaluates 
the required Mellin transforms of the evolution kernels analytically in perturbative form up to 
an order in (aL) k at which sufficient numerical stability is achieved. In the past the analytic 
resummation of the soft- and virtual leading order contributions was performed [9] . In the flavor 
non-singlet case analytic results were obtained to 5th order in aL [10,11] 3 . We compactify this 
description using the small-x resummations of terms oc aln(x), as carried out in [13] to 3rd 
order, to 5th order. However, these terms do not yield the complete description yet. The flavor 
singlet contributions have to be calculated to the same order, which is one of the objectives of 
the present paper. 

The second class of universal corrections concerns the resummation of contributions of the 
type (a\ri 2 (x)) k in the cross section. These terms are resummed by infrared evolution equations. 
The flavor non-singlet contributions were studied earlier in [14,15]. As the comparison of the 
first two terms of the resummation with QED corrections [16] shows, the resummation makes a 
correct prediction for the terms which occur in initial state radiation (ISR), whereas it remains 
an open problem, whether this is the case for final state radiation (FSR) 4 . For this reason we 
apply the resummation only for initial state radiation. In extension of earlier work in which the 
non-singlet resummation was dealt with [15] the resummation is carried out also for the singlet 
contributions. 

The paper is organized as follows. In section 2 we revisit the non-singlet solution to O [(aL) 5 ] 
and compactify previous results w.r.t. exponentiations for small and large values of x. Section 3 
contains the general formalism to derive the singlet solution both for the leading order O (aL) 



terms and those due to the resummation of contributions oc O 



{a\n 2 (x)) h 



consistently. The 



iterated singlet evolution kernels are presented to O [(aL) 5 ] in section 4. In section 5 we derive 

2 Numerical studies in the case of QCD were performed in [7]. For a review see [8]. 
3 Earlier results up to 3rd order were given in [12]. 

4 Note that these terms do potentially violate the Gribov-Lipatov relation, see e.g. [17,18]. 
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the coefficients for the resummation of the O (a\n 2 (x)) k terms for singlet evolution. In sec- 
tion 6 numerical results are presented and section 7 contains the conclusions. Technical aspects 
of the calculation being of interest also for other higher order QED and QCD calculations are 
summarized in the Appendices where convolution tables, which were used in the present calcu- 
lation, are given in Appendix A. As well we summarize a series of higher order relations between 
Mellin transforms and harmonic sums, or polynomials of them, which were used to calculate the 
more complicated convolutions in Appendix B and C. 



2 The Non-Singlet Solution 

The universal radiative corrections discussed in the present paper can be derived using the 
renormalization group equations for mass factorization. We consider a one-flavor system out of 
electrons (positrons) and photons. The flavor non-singlet and singlet distributions are defined 
by the following QED-distributions 



D^ s (a,x) = D e (a,x) - D e+ (a,x) (1) 

D^(a,x) = D e ~(a,x) + D e+ (a,x) . (2) 

Here x denotes the momentum fraction of the collinear daughter-particle radiated off the source, 
< x < 1, and a = a(Q 2 ) = e 2 / (47r) 2 is the running coupling constant with e the electric charge. 
For this system only one non-singlet 'minus' distributions occurs. 

The non-singlet contribution to the polarized QED evolution kernel is the same as in the 
unpolarized case since the splitting functions 

^(*) = 2 (t=t) + ( 3 ) 
are identical. The +- operation is defined by 

£ dx [F(x)] + 4>{x) = £ dxF(x) [<f>(x) - 0(1)] (4) 

for the functions <f>(x), F(x) with x e [0, 1]. 

The leading order non-singlet evolution equation reads 

D e m (a,x) = a(Q 2 )Pg(x) ® D e m (a,x) . (5) 



<91n(Q 2 , 

Here ® denotes the Mellin-convolution 

A(x) <g> B(x) = [ dxi [ dx 2 5(x - x 1 x 2 )A(x 1 )B(x 2 ) . (6) 
Jo Jo 

One may solve (5) using the Mellin-transform 

M[F(x)](N) = f 1 dxx N - 1 F(x) . (7) 
Jo 

The Mellin-transform of the splitting function is the anomalous dimension (N). 
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The QED running coupling constant is given by 



with the first coefficients of the /3-function (3 = —4/3, /?i = —4. The solution of Eq. (8) in 
leading order reads 

<Q 2 ) = — ^- , (9) 
1 - -a L 

with ao = a(^e) an d ^ = l n (*3 2 / m e)- F° r a o-^ ^ 1 & U later results obtained in a may be 
expanded in the parameter a^L to a finite order. 

In Mellin-space the solution of the non-singlet evolution equation reads 



1 + (3 a In — 



g > . (10) 



e / J 



Here we used the boundary condition D{/ s (N,m^) = 1 which corresponds to 

Dg(x,m 2 e ) = 5(l-x) . (11) 

For Eq. (10) one often defines [3 = (2a /7i") ln(<3 2 / m e)- 

Unlike the case of QCD evolution of extended parton densities [19] or structure functions [18], 
which fall rapidly as x — > 1, the input distributions in QED are of the form (11) and are 
distribution-valued at x — 1. Therefore, the Mellin-inversion of (10) cannot be performed by 
a numerical contour integral but has to be worked out calculating all the convolutions up to a 
finite order in x-space in to (aL) k at which sufficient convergence of the series is being obtained. 
In the present paper the iteration is performed up to k = 5. The corresponding convolutions 
are evaluated in explicit form in terms of Nielsen integrals in Appendix A. The calculation does 
partly require the use of multiple harmonic sums [20-24], which can be expressed as Mellin 
transforms of Nielsen integrals. In a previous paper [21] a series of these Mellin transforms has 
been calculated already. In Appendix B additional representations, which were derived and used 
in the present calculation, are given. 

We compared our results to those given previously in [10, 11]. We agree with [10] but find 
a deviation from [11] in the 0((aL) 5 ) term. The result of [11] is very lengthy and can be 
compactified by exponentiating the soft- and virtual contributions to all orders in (aL) k , which 
was done before in Ref. [10]. Moreover, also the latter representation can be shortened in 
exponentiating also the contributions oc (f3\n(x))\ well-known from QCD [25], which has been 
carried out for the first three orders in [13], but is not contained in [10, 11]. 

The non-singlet splitting kernel (no running coupling) can be represented exponentiating the 
soft and small x contributions to all orders and treating the other contributions to finite order 
as 



D NS (x,P) 



exp[(/3/2)(3/4 - lE )\ /3 q _ ^Ii (MM*)) 




:i2) 



r(l + /3/2) 2^ i-p\n(x)Y/i 
with h(z) the associated Bessel function, ^ E the Euler-Mascheroni constant, and 
^o(x) = l+x 2 (13) 
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#i(a;) 
V 3 (x) 



[1 — x) 2 + x 2 ln(x) 



= -(1 - x) [1 - x - x ln(x) + (1 + x)Li 2 (l - x)] 

= "^j 6 ! 1 " ^ 2 )[2Li 3 (l - x) + Li 2 (l - x)} + 5(x - 1) 



(14) 
(15) 
(16) 



+ (1 + 7x 2 )[ln(x)Li 2 (l - x) + 2Si, 2 (l - x)] - Q + 6x - yx 2 ) ln(x) + -^x 2 ln 3 (x) j 



96 



l-x 2 ) 



24Li 4 (l - x) + 12Li 3 (l - x) - -Si, 3 (l - x) 



- 12S 2j2 (l - x) - - ln(x)S 1)2 (l - x) - ^ ln 2 (x)Li 2 (l - x) + 7Li 2 (l - x) 



+ 4(1 + x 2 )Li 2 (l - x) + (1 - 8x + 7x 2 



ln(x)Li 2 (l - x) + 2Si )2 (l - x) 



+ 2(1 + 7x 2 ) ln(x)Li 3 (l - x) - - + 5x - — x 1 ln(x) 



7 1, 2 



23 



I i 

— x(l - x) ln 3 (x) - —x 2 ln 4 (x) + (1 - x) 2 
12 48 



2 + 8 lnW 



(17) 



Except the Oth term all the functions ^k{ x ) behave at least as oc (1 — x) for x — > 1. Therefore in 
these orders the contribution oc 1/(1 — x) in the function D NS (x,/3) is canceled. As mentioned 
before the leading small x terms oc (/31n(x)) fc are resummed into the function Ii (— (3\n(x)) 1 / 2 . 
Indeed many more logarithmic contributions are absorbed in this way from the expanded ex- 
pression. To 3rd order in f3 the functions ^(x) are finite as x — > 0. In the higher orders some 
terms of the order f3 k \n h ~\x), I > remain, but contain a small weight factor only. 



3 The Solution of Singlet Evolution Equations 

In this paragraph we derive the explicit solution of the singlet evolution equation in terms of the 
running coupling constant a = a(Q 2 ), with a = a/(47r) and account both for the leading order 
terms as well the leading terms at small x oc a k hi 2<yk ~ l \x). Since the latter terms stem from 
the higher order evolution kernels, which do not commute with the (resummed) leading order 
solution in general, the perturbative solution has to be constructed accordingly. A framework 
of this kind was derived for Mellin iV-space in Ref. [19] and for the lowest three orders in [26]. 
Here we will work in x-space and have therefore to extend a prescription given to second order 
in [27] to all orders. 

We substitute the scale dependence of the (singlet) evolution equation 5 in terms of a, Eq. (8). 
One obtains 



dD s (a,x) 
da 



1 P (x) + a P^x) + a 2 P 2 (x) + . . . 
a /3 + a Pi + a 2 p 2 ■ ■ ■ 

Pi 



D s {a,x) 



P a 



P (x) + a I Pi(x) - g P (x) + 0(a 2 



Dc(a,x) 



5 The corresponding structures for the non-singlet evolution equations are those of entry (1,1) and are not 
given separately. We will use this structures also resumming the small x terms oc ln 2 (x) later on. 
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1 

a 

The functions Rk( x ) are given by 



RoW + Y, akR k( x ) 
k=i 



D s (a,x). (18) 



Ro(x) = ^Po(x) (19) 
Po 

R k {x) = ^P k (x)-j2w R k-i( x )- ( 2 °) 
Po i= i Po 

The leading order solution is 

D s ,o(a,x) = [exp(- Ro(x)\n(a/a Q ))^\<S) D s (a ,x) = E (a,a ,x) <S> D s (a ,x) , (21) 
where E is the leading order singlet evolution operator. We use the short-hand notation 

[/(« = E^®ti9W, (22) 

with <S>i=i the fc-fold convolution. The distribution at the scale a = do is 

D s (a ,x) = 1 5(1 -x) . (23) 
Higher than leading order contributions are found perturbativly by 

D s (a,x) = U(a,x) E (a,a ,x) <g> U(a Q ,x)~ 1 , (24) 

with 

oo 

U(a,x) = 1 5(1- x) + Y / ak U k (x) . (25) 

k=i 

The matrices U k (x) are obtained recursively from the functions R k (x), 

%(x) = Ri(x) (26) 
_ k-i 

R k (x) = R k (x) + J2 Rk-i® U k (x) . (27) 

i=i 

In general the matrices U k do not commute. The higher than leading order contributions dealt 
with in the present paper concern contributions to the splitting functions oc a k \n 2h (x). If one 
would work in x-space many convolutions of these functions with those representing the fixed 
order solutions .Eo would have to be calculated. The series in a k \n 2h (x) may need more terms 
than in the case of the fixed order solution in 0[(aL) k ]. The specific structure of the terms 
a k \v? k (x) allows to work in iV-space since the Mellin transform 

M[\^ k (x)}(N) = (28) 

is sufficiently damped for k > 1 in the limit Re(iV) — > — oo. Due to this we may calculate the 
t/-matrix in iV-space and use the solution 

D s (a,x) = E (a,a ,x) + U(a,x) <g> E (a,a ,x) + E (a,a ,x) <8> U(a ,x) 

+ U(a, x) <g> E (a, a , x) <g> U(a , x) , (29) 

where 

U(a,x) = U(a,x) -15(1- x) (30) 
U(a,x) = U(a,x)- 1 -15(1- x) . (31) 
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4 The Leading Order Solution to 0[(cxL) 5 ] 

In the present paper we limit the explicit representations to order n = 5 in (a L) 6 . It is 
therefore convenient to expand (21) as 



E (a,x) 



1 5(1 -*) + P (*)a L + 



'f ) (x) + -P (x))M 2 



(32) 



+ ( - r Pf(x) + P^(x) + £ P (x) ) (oqL) 



»(i) 



16 



27 



+ Pf (*) + ^ + | P?>(s) + g P (x)) (a LY 



»(4), 



120 



*) + 5 J*W + ^ i*W + 1 *"?'<*) + H fl><*) ) ( if 



Eq. (32) makes the effect of the running coupling (9) explicit. The matrices P^ are 

P ( k \x) = ®t 1 P (x) . (33) 

The corresponding expressions for the (1, 1) -components of P^\x) are given relative to the 
non-singlet components P^g(*) of section 2. The singlet-matrix components Pjj (*) given 
below (34-53) were calculated using the convolution formulae of Appendix A. 7 The projections 
P^j describe the splitting of an electron into an electron (1,1), of a photon into an electron, 
positron, respectively (1,2), an electron into a photon (2,1), and a photon into a photon (2,2) in 
kth. order in the renormalized coupling constant. 

The individual components for the different orders are given by : 



0th order 



P[%) = -4(1-2*) 
Pff(x) = 2(2-*) 



--5(1 -x). 



(34) 
(35) 
(36) 

(37) 



Here a non-vanishing contribution to P^fa) emerges due to 4-momentum conservation, see 
also [29]. 



1st order : 



pj»( x ) = pW + 8[5(l-*)+2(l + *)ln(*)] 

P§(x) = 4 

P$(x) = 4 

P$(s) = 8 



2(2*- l)[21n(l -*) -ln(*)] + 
(2 -*)[21n(l -*) + ln(*)] - - 



13-8* 
3 

- 13* 



6 



5(1 - *) + 2(1 + *) ln(*) + -5(1 

9 



* 



(38) 
(39) 

(40) 

(41) 



6 Beyond this level we will use only the exponentiation for the non-singlet terms. 

7 Note that the second order contribution P^(x) given in the literature is partly incorrect [28] for polarized 
scattering. 
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P$ + 16 



49 4 
- — (1 - x) + 20(1 - x) ln(l - x) - -(7 - 8x) ln(x) 



8 



+8(1 + x) ln(x) ln(l - x) - 2(1 + x) ln 2 (x) + 8(1 + x)Li 2 (l - x) 

1379 - 1318a; n/i n . . . 44 - 16a: , ,„ , 98 + 1583;, . . 

— + 8(1 - 2x)C(2) + ln(l - x) + ln(x) 



18 ' " /,v ' ' 3 ' 3 

-2/i \ , o/i n«\l„/«\ 1„/1 „\ , o/i n„\,„2/ 



-8(1 - 2a;) ln 2 (l - x) + 8(1 - 2x) ln(x) ln(l - x) + 3(1 - 2x) ln 2 (x) 

659 1379 J/rt 8 -22a;, „ , 71,, x , , N 
-— + — x - 4(2 - x)C(2) — In(l - x) - — (1 + x) ln(x) 

3 

+4(2 - x) ln 2 (l - x) - 4(2 - x) ln(x) ln(l - x) - -(2 - x) ln 2 (x) 

79 20 
■ — (1 - x) + 20(1 - x) ln(l - x) - —(2 - x) ln(x) 

o o 



(42) 



(43) 



(44) 



(45) 



+8(1 + x) ln(x) ln(l - x) - 2(1 + x) ln 2 (x) + 8(1 + x)Li 2 (l - x) - —6(1 - x 

— I 



- (*|* + 120C(2)) (1 - x) - ^(1 - x) ln(l - x) 



- ( 1445 + 9 3089X + 48(1 + x)C(2)) ln(x) + 120(1 - x) ln 2 (l - x) 

_^!(4 _ 5a; ) i n ( x ) i n (i _ x ) _ °° l n 2 (x) - 24(1 + x) ln 2 (x) ln(l - x) (46) 

3 3 

+48(1 + x) ln(x) ln 2 (l - x) - ^(1 + x) ln 3 (x) + ^(1 + x)Li 2 (l - x) 
+96(1 + x) ln(l - x)Li 2 (l - x) - 96(1 + x)Li 3 (l - x) + 48(1 + x)Si, 2 (l - x) 



16 



42587 10768 124 - 32rr > .„, nn/ , n _, oN 
+ x o C(2) - 32(1 - 2x)C(3) 



108 27 

/2693- 2506x . x>/ A . 
+ f g + 48(1 - 2x)C(2) J ln(l - x) 

/3389 - 1384x A/ x>/ A . , 124 - 32x , 2/ 
- f + 24(1 - 2x)C(2) j ln(x) + ln 2 (l - x) (47) 

380 + 6UX ln(x) ln(l - x) - m + 3MX ln 2 (x) - 16(1 - 2x) ln»(l - x) 



3 v ' v ' 6 

+24(1 - 2x) ln(x) ln 2 (l - x) + 10(1 - 2x) ln 2 (x) ln(l - x) 

--(1 - 2x) ln 3 (x) + (168 + 204x)Li 2 (l - x) + 44(1 - 2x) ln(x)Li 2 (l - x 
o 

+52(l-2x)S lj2 (l-x) 
_ T5384 42587 16 - 62x ,„, rt 

16 VW + ~zm x + — 3 — c(2) + 16(2 " x)c{3) 

(12^3 2693 \ 
— — x + 24(2-x)C(2)jln(l-x) (48) 



/907 1997 v/n A, . 

+ 1 — —x + 12(2 - x)((2) ) ln( 



16 - 62a; 



9 36 
290 + 314a; 



ln 2 (l -x) 



ln(x) ln(l — x) + 



y + ln2 (^) + 8 ( 2 - x ) ln3 (! " 

-12(2 - x) ln(x) ln 2 (l - x) - 5(2 - x) \n 2 (x) ln(l - a;) + -(2 - a;) ln 3 (a;) 

6 

-22(2 - x) ln(x)Li 2 (l - x) - (102 + 84x)Li 2 (l - x) - 26(2 - x)Si, 2 (l - x) 
/2177 V , 236 



16 



{ , +40C(2)J(l-a;)-— (1 - x) ln(l - x) 



785 



- ( 183 + ^x + 16(1 + x)((2) J ln(ar) + 40(1 - x) ln 2 (l - x) 



(49) 



128 - 112a; 



ln(x) ln(l - x) 



103 - 107a; 



ln 2 (x) 



+ 16(1 + x) ln(x) ln 2 (l - x) - 8(1 + x) \n 2 {x) ln(l - a;) - 2(1 + x) ln 3 (x) 

o 

+32(1 + x) ln(l - x)Li 2 (l -x)- -(1 + x)Li2(l - a;) - 32(1 + x)Li 3 (l - a;) 
+16(1 + x)Si )2 (l-x) + ^5(1 -x) 



4th order 



P^g + 32 
/15212 



/ 73678 



27 



I- 



+ 704C(2) + 640C(3) (1 - x 



9 



+ 960C(2) (1 - x) ln(l - x) 



+ 



+ 



^36568 90C8 



27 27 
704(1 -x) ln 2 (l -x) 



x + 416C(2) - 544((2)x + 256(1 + x)((3) ) ln(x) 
/8000 + 20672a; 



9 



+ 384(1 + x)C(2) ln(x) ln(l - a;) 



/ 2696 + 4472a; 
V 9 



+ 96(1 + x)) \n 2 (x) + 320(1 - x) ln 3 (l - x) 



-(416 - 544x) ln(x) ln 2 (l - x) - (152 - 88a;) In 2 (a;) ln(l - x) 
40 

+ — (8 - 7x) \n 3 {x) + 128(1 + x) ln(x) ln 3 (l - x) - 96(1 + x) ln 2 (a;) ln 2 (l - x) 



(50) 



9 

+ \(l+x)\n\x) 



f 14336 
V 9 



+ 384C(2)) (l + x)Li 2 (l-x) 



+128(1 + x) ln(l - x)Li 2 (l -x)- 720(1 - x) ln(x)Li 2 (l - x) 

+384(1 + x) ln 2 (l - x)Li 2 (l - x) - 96(1 + x) ln 2 (x)Li 2 (l - x) 

-128(1 + x)Li 3 (l - x) - 768(1 + x) ln(l - x)Li 3 (l - x) + 768(1 + x)Li 4 (l - x) 

-(816 - 944x)Si )2 (l - x) + 384(1 + x) ln(l - x)S li2 (l - x) 

-256(1 + x) ln(x)Si, 2 (l - x) - 384(1 + x)S 2)2 (l - x) - 192(1 + x)Si, 3 (l - x) 
32 



1912189 958405 7852 - 7064a; >/ . 640 - 128x A/x t/ x>/ , x 
x : C(2) + - C(3) - 48(1 - 2x)C(4) 



648 



324 



9 



9 



+ (-^r + ^r x ~ 64(5 - xK(2) - 256(1 - 2;r)c(3) ) ln(1 - x) 

+ - - 368C(2) - 656C(2)x + 128(1 - 2x)C(3)) ln(x) 

+ ( 7852 " 9 7064 % 192(1 -2,) C (2)) 1^(1-,) (51) 

- ( 922 °~ 4616X + 192(1 - 2x)C(2)) ln(x) ln(l - x) 

+ (-^ + ^ - 24(1 - 2s)C(2)) ln 2 (x) + |(5 - x) ln 3 (l - x) 



+ (368 + 656x) ln(x) ln 2 (l - x) - 616 + 1000x ln 2( x ) ln(1 _ x ) 

o 

- (^20 + yx) ln 3 (x) - 32(1 - 2x) ln 4 (l - x) + 64(1 - 2x) ln(x) ln 3 (l - x) 

+24(1 - 2x) ln 2 (x) ln 2 (l - x) - — (1 - 2x) ln 3 (x) ln(l - x) - —(1 - 2x) ln 4 (x) 

3 X 2 

-(152 + 272x)Li 2 (l - x) + (1056 + 1248x) ln(l - x)Li 2 (l - x) 
32 

- — (4 + x) ln(x)Li 2 (l - x) + 288(1 - 2x) ln(x) ln(l - x)Li 2 (l - x) 

-16(1 - 2x) ln 2 (x)Li 2 (l - x) - (1056 + 1248x)Li 3 (l - x) 

-288(1 - 2x) ln(x)Li 3 (l - x) + 352(1 - 2x) ln(l - x)Si, 2 (l - x) 

+80(1 - 2x) ln(x)Si, 2 (l - x) - 32(1 - 2x)Li 2 (l - x) + 48(9 + 13x)S li2 (l - x) 



-224(1 - 2x)S 2i2 (l - x) + 224(1 - 2x)S 1>3 (l - x 
32 



958405 1912189 3532 - 3926x 64 r x>/oX ni . n 

x + C(2) - — (1 - 5x)C(3) + 24(2 - x)C(4) 



648 1296 9 3 

+ ( 41849 + 4 94913x + 296C(2) + 344C(2)x - 64(2 - x)C(3) ) ln(x) 



/29182 - 28640x . „ , , nn/n , //0 A, .„ 

+ [ + 32(1 - 5x)C(2) + 128(2 - x)C(3)J ln(l - x) 

/ 3532 - 3926x nn . n .\ , 

+ ( g 96(2-x)C(2)jln 2 (l-x) (52) 

+ ~ 9 3M2X + 96(2 - x)C(2)) ln(x) ln(l - x) 

+ - 6 ~w x + 12(2 - x)c(2) ) ln2(x) - f (1 - 5x) ln3(1 - x) 



, . 436 + 484x, 9 . . , ,„ . /35 14 * , o. , 

-(296 + 344x) ln(x) ln 2 (l - x) + ln 2 (x) ln(l - x) + I — + — x ) ln 3 (x) 

3 V 2 3 / 

+16(2 - x) ln 4 (l - x) - 32(2 - x) ln(x) ln 3 (l - x) - 12(2 - x) ln 2 (x) ln 2 (l - x) 

20 5 
+ — (2 - x) ln 3 (x) ln(l - x) + —(2 - x) ln 4 (x) + (136 + 76x)Li 2 (l - x) 

16 

-(624 + 528x) ln(l - x)Li 2 (l - x) (1 + Ax) ln(x)Li 2 (l - x) 

3 

-144(2 - x) ln(x) ln(l - x)Li 2 (l - x) + 8(2 - x) ln 2 (x)Li 2 (l - x) 

+48(13 + llx)Li 3 (l - x) + 144(2 - x) ln(x)Li 3 (l - x) + 16(2 - x)Li 2 (l - x) 

-312(1 + x)Si, 2 (l - x) - 176(2 - x) ln(l - x)Si, 2 (l - x) 

10 



-40(2 - x) ln(x)Si, 2 (l -x) + 112(2 - x)S 2 , 2 (l - x) - 112(2 - x)Si, 3 (l - x 
P$(x) = 32 f^^ + ^C(2) + 160C(3)Vl-x) -(1417 + 240C(2))(l-x)ln(l-x) 



108 3 
/ 32002 - 6257a; 352 - 368a; , 

+ ( 27 + 3 C(2) + 64(1 + x)C(3) ' lH( ' r ] 

628,, N1 „, N /2150 + 3406x 



^(1 - x) ln 2 (l - x) - ( + 3 + 96(1 + x)C(2)) ln(x) ln(l - x) 



+ ^ 763 + 626x + 24(1 + x)C(2)) ln 2 (x) + 80(1 - x) ln 3 (l - x) (53) 
352 " 368a l ln(x) ln 2 (l - x) - 394 - 38to ln 2 (x) ln(l - x) 



3 
50 

+ — (5 - 4x) ln 3 (x) + 32(1 + x) ln(x) ln 3 (l - x) 
9 

20 7 
-24(1 + x) ln 2 (x) ln 2 (l - x) - — (1 + x) ln 3 (x) ln(l - x) + -(1 + x) ln 4 (x) 

3 6 

16 

-(926 + 96C(2))(1 + x)Li 2 (l - x) + — (1 + x) ln(l - x)Li 2 (l - x) 

-380(1 - x) ln(x)Li 2 (l - x) + 96(1 + x) ln 2 (l - x)Li 2 (l - x) 

16 

-44(1 + x) ln 2 (x)Li 2 (l - x) - — (1 + x)Li 3 (l - x) 

-192(1 + x) ln(l - x)Li 3 (l - x) + 192(1 + x)Li 4 (l - x) 
1 9^2 — 1 268r 

S 1>2 (1 - x) + 96(1 + x) ln(l - x)S 1>2 (l - x) 

-104(1 + x) ln(x)Si )2 (l - x) - 96(1 + x)S 2)2 (l - x) - 88(1 + x)Si, 3 (l - x) 

:!2 -5(l-x) 
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We expressed the above splitting functions in terms of poly-logarithms and Nielsen integrals [30] , 

/_1 \n-l „i f j 7 

Li n (x) = 5 n _ M (x) = / -ln«- 2 (^)ln(l-x^) (54) 

[n — Z)\ Jo z 



(n-2)!io z 

i 

(n- l)\p~ 

5 Oi n+1 ln 2n (cc) Corrections to The Evolution Kernel 



S n , p (x) = ^lyln^im-^). (55) 



The resummation of the leading small-x terms to all orders in the coupling constant were con- 
sidered in QCD in Refs. [7, 14] and [6] for the singlet case. Very recently the results of [5] in 
the interpretation of Ref. [6a, 12] were verified to be correct in the 3rd order in the coupling 
constant in the unpolarized non-singlet case [32]. This gives further confidence in the method 
used. The QED analogue is obtained by setting Cp — T F — 1,Ca — 0. We will furthermore 
consider only one flavor, the electron. The corresponding contribution to the matrix of singlet 
splitting functions reads 

OO 1 

P(x,a)^ = £ P iV +1 m 2i (a;) = ~^M~ X [F (N,a)\ (x). (56) 
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The matrix F (N, a) obeys the equation 

a , , 8a 



Fo(X,,) = | () ;r 2 ^M -^F 8 (iV,a)G + ^iF2(iV,a), 



where 



9 a , 1 1 ., 



F 8 (iV, a) = 16^-M 8 + ^^(iV, a). 



The generating matrices are 
Since Mg = M 8 in the case of QED, F 8 (iV, a) can be given in analytic form 



1 "M 

o or 



Go 



1 




z 


pl'J 


^fg 


p(0 
gg 





0.2000000000E+01 


-0.4000000000E+01 


0.0000000000E+00 


1 


-0.1400000000E+01 


-0.4000000000E+01 


-0.8000000000E+01 


2 


-0.8666666667E+01 


0.6666666667E+01 


-0.2666666667E+01 


3 


0.1444444444E+01 


0.2444444444E+01 


0.2044444444E+01 


4 


0.7888888889E+00 


-0.2825396825E+00 


0.4634920635E+00 


5 


0.5537918871E-02 


-0.1008112875E+00 


-0.3626102293E-01 


6 


-0.1216503661E-01 


-0.5536849981E-03 


-0.9307893752E-02 


7 


-0.5063924905E-03 


0.8396483000E-03 


-0.3947845218E-04 


8 


0.4517396184E-04 


0.2736857499E-04 


0.4555056936E-04 


9 


0.2611470148E-05 


-9.1960462641E-05 


0.1192943787E-05 


10 


-0.4136771391E-07 


-0.9306140011E-07 


-0.6994805938E-07 


11 


-0.4259181787E-08 


0.1228221436E-08 


-0.2774664782E-08 


12 


-0.1716660622E-10 


0.1073980741E-09 


0.3096908050E-10 


13 


0.2808782930E-11 


0.3728781778E-12 


0.2327857737E-11 


14 


0.3713452748E-13 


-0.5274589147E-13 


0.7034915420E-14 



Table 1: The coefficients of the matrices P^ ' n Eq. (56). 
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Note that Eq. (60) has no pole singularities but a branch cut for complex values of N. Due to 
this only moderate effects are expected. The coefficient matrices P^o are obtained by iterative 
solution of Eqs. (57). The first coefficients read 

P™„ = l(j "J) (61) 
Piio = 2 ( ~\ Zl ) (62) 

= I ( :S -4° ) • w> 

These matrices agree with the constant contribution to the leading order polarized singlet split- 
ting function and the coefficients of the In 2 (a;) contribution to the next-to- leading order splitting 
function [31] adjusting the color factors as above. As already noted in the case of QCD the 
off-diagonal matrix elements are related by 

Pg = -P?}. (64) 

The first 15 coefficients are listed in Table 1. The corrections oc a n+1 \n 2n (x) are due to the 
most singular parts of the non-leading order anomalous dimensions and are universal because 
they emerge as part of a conformal solution. For each power in a we have to construct the 
corresponding matrices C7 fc , Eq. (25). We use Eqs. (30, 31) and work in iV-space. The Mellin- 
inversion is carried out by a numerical contour integral around the singularities in the complex 
iV-plane. Finally the solution is given by Eq. (29). 



6 Numerical Results 

The calculation of the universal radiators is carried out to a precision which will not need further 
improvement in the range of collider energies at present facilities and those in the foreseeable 
future. These corrections need to be supplemented by the non-universal contributions of the 
respective processes under consideration. In various cases the universal contributions dominate 
in some kinematic regions. They do, however, not substitute any complete calculation. In the 
above we presented the respective expressions in analytic form which are easily implemented 
into QED radiators used in analysis programs. Now we turn to numerical illustrations. 

The universal radiator function D NS (x, Q 2 ), Eq. (12), summing the leading log corrections up 
to 0((aL) 5 ) is shown in Figure 1 as a function of the momentum fraction x for different values 
of Q = ^fQ 1 = 10, 100, 1000 GeV, the range of relevant scales at high energy colliders. The 
function rises strongly for growing values of x and reaches 1 for x ~ 0.9. D^six, Q 2 ) grows with 
Q, however, the scale dependence is weaker. In Figure 2 we compare the soft photon resummation 
beyond the 5th order with D^ s (x, Q 2 ). Although these terms become more significant for large 
values of x the correction is widely below 3 ppm in the range of Q considered above, which 
shows that summing the leading orders to 0((aL) 5 ) for the full radiator is already a sufficiently 
accurate approximation. In Figure 3 we show the relative impact of the radiator in first order 
compared to the resummed radiator for the same values of Q. This ratio varies between 90 
to 95 % for small values of x to 120 to 130 % for x — > 1. Also the the 4th and 5th order 
approximation for D k ' ies (x,Q 2 ) are compared, including the respective soft exponentiation. The 



13 



deviation is smaller than 1.5 ■ 10~ 5 at small values of x and vanishes for large values of x, showing 
the degree of convergence. 

In Figure 4 we show the singlet contributions for the leading logarithmic radiators in case of 
longitudinally polarized electrons or positrons. Here, Dn(x,Q 2 ) denotes the pure singlet part 
for the fermion-radiator to which the non-singlet contribution has to be added. The spread 
of the functions w.r.t. the scale variation in Q is larger for the off-diagonal radiators D i2 and 
D 2 i, as well as for D 22 if compared to the fermionic radiator D\\. The off-diagonal radiators 
Du(x) and D 2 i(x) vary between —6 and +6 %, and +7 and +2 %, respectively from low to 
large values of x. The pure singlet part of the diagonal radiators is smaller. It grows from 
values of —0.2 %(— 0.1 %) for small values of x to +0.05 % and vanish as x — > 1. The relative 
smallness of these contributions results from the fact that they start with 0((aL) 2 ) only aside 
of contributions oc 5(1 — x) not shown in the Figures. 

Figure 5 compares the size of the first order contribution to £) u , D 12 and D 22 for the chosen 
values of Q with the total contribution. For Z? n the first order contribution is 20 % larger than 
the total contribution at small x while it basically yields the full contribution for large values of 
x. The first order contribution to Di 2 dominates Di 2 , except of the small region x ~ 0.5, where 
higher order contributions are significant. D 2 i receives higher order corrections. D^/D^ — 1 
varies from —5 % to +15 % from small to large values of x. In Figure 6 we compare the 
convergence of the approximation showing the impact of the 5th order term w.r.t. the first four 
orders. This effect is of O(10~ 5 ) with some variations, where the effect in the case of D 22 is an 
order stronger for large values of x. This shows that the this level of resummation is sufficient. 

In Figure 7 we compare the radiators due to the resummation of the a ln 2 (x)L-terms for the 
same choice of scales Q as before. D n (x, Q 2 ) denotes the complete fermionic radiator. The non- 
singlet contributions were dealt with before in Ref. [15]. The pure singlet term contributes only 
with 0(a 2 ). Comparing Table 1 with Table 1 [15] one sees that the pure singlet contributions 
change the non-singlet result and is not subleading. We display the contributions starting with 
0(a 2 ), i.e. those containing at least one term oc \n 2 (x). The first order contribution, the small 
x limit of the leading order splitting functions, has been contained in the radiators of the fixed 
order calculation already. The respective part of the radiators \Dij(x,Q 2 )\ is of 0(1 %) at small 
values of x and vanishes towards large values of x. 

We finally apply the corrections to polarized deeply inelastic charged lepton scattering off 
polarized protons as an example. The hadronic tensor of the polarized part of the eiV-scattering 
cross section in case of pure photon exchange 8 reads 

W/V = ttiw>v-^j-9i(x, Q ) + is^xo ^—y 2 g 2 {x, Q ) . (65) 

S a denotes the nucleon spin vector, p the nucleon momentum, and q the vector of the 4- 
momentum transfer, with Q 2 = —q 2 and x = Q 2 /(2p.q). The polarized part of the scattering 
cross section for longitudinal nucleon polarization Sl, integrated over the azimuthal angle 0, is 



d 2 a(X, ±S L 



B a 2 



= +27TS-- 

dxdy Q 4 



-2Xy [2-y- J x 9l (x, Q 2 ) + 8A^%(*, Q 2 ) 



(66) 



Here M is the nucleon mass, S the cms energy, a the fine structure constant, y = 2p.q/S, A 
is the degree of lepton lepton polarization, x and y are the Bjorken variables, and gi, 2 (x, Q 2 ) 

8 The corresponding expressions in the case of additional weak boson exchange are given in [33]. 
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are the two electromagnetic polarized structure functions, which are represented at twist-2 level 
referring to the parameterization [34]. Very similar results are obtained using other recent 
parameterizations [35,36]. We consider the case of initial state radiation and sum over the final 
state bremsstrahlung. Furthermore, we cut for the Compton peak, [3]. The kinematic variables 
chosen are the variables for leptonic corrections. The rescaled variables x, y, and s are 

x = — -, y = , s = zS , (67) 

z + y — 1 z 

and the radiation threshold is z = (1 — y)/(l — xy) . The radiative correction to the differential 
cross section (66) is given by : 



rl 2 a isr rPn B /■! "> d 2 rr B 

Hi Ju U' U Uj<uU Li' (J J IJ || \Jj<uU Li' (J 



=x,y=y,s=s <* 



(68) 



J(z,y,z) 



(69) 



-D(z, Q 2 ) denotes the respective radiator function for fermion-fermion transitions for z < 1, and 
the Jacobian J is given by 

dx/dx dy/dx 
dx/dy dy/dy 

Previously, the leptonic radiative corrections for this process have been calculated in 0(a) com- 
pletely in Ref. [37], where the complete calculation was also compared to the leading log result 
in 0(aL). In the case of longitudinal polarization studied in the present paper, both at HER- 
MES energies, = 7.4 GeV, and for ep collider energies, = 314 GeV, the complete 0(a) 
corrections deviate at most by a few % for small values of x and very large values of y from the 
leading log result. 

In Figure 8 the initial state QED corrections for the differential scattering cross section (66) 
are shown for the kinematics of the HERMES experiment at DESY with v^S* = 7.4 GeV, as a 
ratio to the double differential Born cross section using the resummed radiator as a function 
of y for characteristic values of x. The radiative corrections grow with y and towards smaller 
values of x and may reach values of O(80 — 100 %) at large y. The effect of the higher than 
leading log corrections to those of leading order, normalized to the differential Born cross section 
is illustrated in the right figure. Depending on the value of x, the corrections in the lower y 
range are of up to O(±0.5 %) and grow rapidly in the high y region reaching 0(4 %). 

For a hypothetical future ep collider operating both with longitudinally polarized electrons 
and protons at a cms energy of v^S* = 1 TeV the double differential radiative corrections due 
to initial state radiation are shown in Figure 9. For fixed values of x all the corrections grow 
strongly towards large value of y. For large values of x they are negative for smaller values of 
y and are everywhere positive for x values below x ~ 0.1 . Due to the larger kinematic range 
the scaling violations due to the running fine-structure constant, a(Q 2 ), are larger than in the 
case of the HERMES experiment. Even for x values of O(10 -2 ) radiative corrections of +160 % 
may be reached for y ~ 0.9. We compare the impact of the higher than leading log corrections 
to those of leading order, normalized to the differential Born cross section. Depending on the 
value of x, the corrections in the lower y range are of up to 0(±2.5 %) and grow rapidly in the 
high y region reaching O(20 %). 

7 Conclusions 

Mass factoriztion at the one hand and infrared evolution equations on the other hand, allow 
to resum two different classes of universal QED corrections oc (a\n(s/m 2 e ) k and a k+l \n 2k (z), 
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respectively, for characteristic invariant masses s and radiation fractions z. The corresponding 
radiator functions are single-particle quantities and describe the universal part of the transitions 
between the light species, electrons e and photons 7, via amplitudes , (i, j) = e, 7, for collinear 
kinematics, similar to the parton model in QCD. In this way these radiators resum the universal 
part of the QED parton cascades. We calculkated the non-singlet and singlet contributions for 
the case of polarized scattering to 0((aL) 5 ) and resuummed the contributions oc a k+l \n 2h (z). 
Both radiators are calculated to an accuracy which is sufficient for the foreseeable applications 
in present and future high energy experiments. While the calculation of the resummation of the 
a k+1 ln 2fc (z)-terms can be done using Mellin space techniques, since the functions contributing 
are damped towards large Mellin parameters Re(iV) — > —00, this is not the case of for the 
radiator, resumming the (a ln(s/mg) fc -terms. Due to this the corresponding convolutions had 
to be calculated analytically in terms of Nielsen functions. The a k+1 ln 2fc (z)-terms are leading 
small z pieces of the higher order splitting functions. To deal with them in the context of the 
resummation of the (a ln(s/mg) fc -terms, one has to invoke the mechanism of general higher order 
solutions of the singlet QED evolution equations for fermionic and photonic densites. Since the 
respective matrices of the splitting functions of the different orders in the coupling constant do 
not commute, it is practically excluded to write the solution in terms of analytic functions in 
z-space. However, one may solve the problem numerically and adds these contributions to the 
analytic radiators found for the (a ln(s/mg) fe -terms. The main objective of the present paper 
was to provide all necessary radiators and to simplify results worked out earlier. The radiators 
obtained can easily be adopted for experimental analysis and simulation programs. 
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8 Appendix A: Mellin Convolutions: Basic Functions 



In this appendix we list the convolutions of functions up to depth 5 9 , which were calculated 
recursively in explicit form, since they are of general interest for other higher order calculations 
in QED and QCD. Some of the integrals require to use Mellin transforms and algebraic relations 
between the finite harmonic sums [21], being associated to them, to be solved. The respective 
expressions which were used in this calculation are given in appendix B and C, as long they were 
not contained in Ref. [21]. 



f(x)®6(l-x) 
1® 1 

1 <g) x 
x <S> x 
1 



XJ + 



X 



1 



1 — X 

1 

1 — X 



'ln(l — x) 



1-x 



1-x 



dn(l — x) 
^x~ 

ln(l — x) 
1 — x 



x 



(rb) + 0ln(1 -^ 

(t-—) ®xh\(l-x) 
VI - x/ + 



X 



1 ' <g> ln(x) 
x ln(x) 

ln(x) 
1 — x 
1 <g) ln(l - x) 
1 <g) x ln(l — x) 
x ® ln(l — x) 
x ® x ln(l — x) 



1-x 
1 

1 — X 



f(x), V f(x) e S'[0, 1] 

-ln(x) 

1-x 

— xln(x) 

ln(l — x) ln(x) 



x 



+ 



x = 



1-x 1 
ln(l — x) — In (a;) 

x [ln(l - x) - ln(x)] + 1 - x 



31n 2 (l -x) 
2 1-x 



1-x 



1-x 



Li 2 (x)-C(2) + -ln 2 (l-x) 



x 



Li 2 (x) - C(2) + - ln 2 (l - x) - ln(l - x) + ln(x) 
2 



+ ln(l -x) 
ln 2 (l - x) - ln(x) ln(l - x) - £(2) 



x 



ln 2 (l - x) - ln(x) ln(l - x) - C(2) 



+(1 -x)[m(l -x) - 1] 

-- ln 2 (x) + Li 2 (l - x) + ln(x) ln(l - x) 

x — - ln 2 (x) + Li 2 (l — x) + ln(x) ln(l — x) 

— 1 + x — x ln(x) 

, , x ln(l-x) lln 2 (x) 

ln(x) — 

w 1-x 2 1-x 

Li 2 (x)-C(2) 

(1 - x) [ln(l - x) - 1] 

xln(x) + (1 — x) ln(l — x) 

x[Li 2 (x)-C(2)] 



(70) 
(71) 
(72) 
(73) 

(74) 

(75) 
(76) 
(77) 
(78) 



(79) 
(80) 

(81) 
(82) 

(83) 
(84) 

(85) 
(86) 
(87) 



3 A series of other useful convolutions has been given in Ref. [18] recently. 
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1 ln(x) 
1 — X 

ln(x) 

x® 

1 — X 

1 <S> ln(x 

1 <S> x ln(x 
x <S> ln(x 

x ® x ln(x 

1 <g> ln 2 (x 

x (g> In 2 (a; 
1 (g) x In 2 (a; 

x ® x ln 2 (x 

1 <g> ln 3 (x 

x (g) ln 3 (x 
1 <g) x ln 3 (x 

x ® x ln 3 (x 

1 <g> ln 4 (x 

x ® ln 4 (x 
1 <g> xln 4 (x 

x ® x ln 4 (x 
ln(x) ln(l — x 



-Li 2 (l-x)--ln 2 (x) 



x ln(x) ln(l — x) 



= x 



-Li 2 (l-x)--ln 2 (x) 



+ 1 — x + ln(x) 



1 



-^ln 2 (x) 

— 1 + x — xln(x) 
1 — x + ln(x) 

-ixln 2 (x) 

-iln 3 (x) 
3 V ; 

ln 2 (x) + 2(ln(x) - x + 1) 

—x ln 2 (x) + 2(x ln(x) — x + 1) 

-ixln 3 (x) 
3 v 1 

--\n\x) 
4 

ln 3 (x) + 3 ln 2 (x) + 6(ln(x) - x + 1) 

— xln 3 (x) + 3xln 2 (x) — 6(xln(x) — x + 1) 

— -xln 4 (x) 
4 



(89) 

(90) 

(91) 

(92) 
(93) 

(94) 

(95) 

(96) 
(97) 

(98) 

(99) 

(100) 
(101) 

(102) 
(103) 



'ln 2 (l - x) 
1 — x 



-iln 5 (x) 

5 

ln 4 (x) + 4 ln 3 (x) + 12 ln 2 (x) + 24(ln(x) - x + 1) (104) 

— xln 4 (x) +4xln 3 (x) — 12xln 2 (x) 

+24(xln(x) -x+ 1) (105) 

-ixln 5 (x) (106) 
5 

2Si, 2 (l - x) - Li 3 (l - x) 

+ [ln(x) + ln(l -x)]Li 2 (l - x) 

~ ln 2 (x) ln(l - x) + ln(x) ln 2 (l - x) 

-C(2)ln(x) (107) 

x|2S li2 (l - x) - Li 3 (l - x) 

+ [ln(x) + ln(l -x)]Li 2 (l - x) 

~ ln 2 (x) ln(l - x) + ln(x) ln 2 (l - x) - C(2) ln(x)| 

+2(1 — x) — x ln(x) ln(l — x) + x ln(x) 

-(1 - x) ln(l - x) - Li 2 (l - x) (108) 

^41n 3 (l-x) , . ln 2 (l-x) 

mix) 

3 1-x V ; 1-x 
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-x/ + 



1 — X 



ln 2 (l-x) 



-^—] <g>xln 2 (l - x) 
— xj + 



1 -x 

1 



® ln 2 (x) 



1 -x 



ln(x) ln(l — x) 
+ 1 — x 



ln 2 (x) 



1 — X J -L- 1 — X 



1 ® ln(x) ln(l — x) 
1 <S> xln(x) ln(l — x) 

1 <g) Li 2 (l - x) 
1 ® xLi 2 (l — x) 

1 ®ln 2 (l -x) 
1 <g>xln 2 (l - x) 



-2C(2)!^ + 2 C (3). ! 



®Li 2 (l-x) = 

1 - x/ + 



® xLi 2 (l — x) = 



J + 



1 — x 1 — X 

-Si j2 (l - x) - Li 3 (l - x) 
+Li 2 (l - x)[ln(l - x) - ln(x)] 
-x [Si, 2 (l - x) + Li 3 (l - x) 

-In (^^) Li 2 (l-x) 
-(1 - x)[l - Li 2 (l - x)] - xln(x) 
ln 3 (l - x) - ln(x) ln 2 (l-x) 
-2C(2)ln(l-x) + 2C(3) 



x 



ln 3 (l - x) - ln(x) ln 2 (l-x) 



-2C(2)ln(l-x)+2C(3) 
+ (1 - x) [ln 2 (l - x) - 2 ln(l - x) + 2 
-2 Li 3 (x)-C(3)-ln(x)C(2) + iln 3 (x) 



xln 2 (x) = -2xLi 3 (x) + 2[C(3)-l]x + 2[C(2) + l]xln(x) 



-xln 2 (x) xln 3 (x) + 2 

3 

Si, 2 (l - a:) ln(x)Li 2 (l - x) _ ln(x) 
1-x 1-x ~ U J l-x 

lln 2 (x)ln(l -x) ln(x)ln 2 (l - x) 



= 2 



2 1-x 
S 1>2 (l-x) lln 3 (x) 



1-x 



1 — x 



3 1 -x 



_ 2 Li 3 (x)-C(3) + ln(x| 

1 — X 1 — X 

C(3) + Li 2 (x) ln(x) - Li 3 (x) 

-Li 2 (l - x) + (1 - x)[2 - ln(l - x)] 

+xln(x)[l -ln(l -x)] 

-2[S li2 (x)-C(3)] + iln 3 (l-x) 

-2Si, 2 (l - x) - ln(x)Li 2 (l - x) 
(1 — x)Li 2 (l — x) — 1+x — x ln(x) 
2[C(3) - S 1)2 (x)] 

(1 - x) [ln 2 (l - x) + 2 - 2 ln(l - x) 
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1 In 2 (a;) 

1 — x 
ln(x) ln(l — x) 
1 -x 



x <S> ln(x) ln(l — x) 

x <S> x ln(x) ln(l — x) 
ln 2 (l -x) 



x <g> 



1 -x 



x <S> 



x <g> Li 2 (l — x) 

x <g> xLi 2 (l — x) 
x <S> ln 2 (l - x) 
x (g) x ln 2 (l — x) 
ln(x) ln(l -x) 



x 



= 2S lj2 (l-x)-iln 3 (x) (125) 
= Si )2 (l - x) - Si j2 (x) 

+ - ln(x) ln(l - x) [ln(l - x) - ln(x)] + C(3) (126) 

= xLi 2 (l -x) + -ln 2 (x) + (1 -x)ln(l -x) 

+ ln(x)[x + ln(l -x)] (127) 

= x[C(3) + ln(x)Li 2 (x) -Li 3 (x)] (128) 

= 2x[C(3)-S li2 (x)] + (l-x)ln 2 (l-x) 

+2x[C(2)-Li 2 (x)] + ^ln 3 (l-x) (129) 

= (1 -x)Li 2 (l -x) - ixln 2 (x) (130) 

= -x[2S 1)2 (l-x)+ln(x)Li 2 (l-x)] (131) 

= (1 -x)ln 2 (l -x) -2x[Li 2 (x) -C(2)] (132) 

= 2x[C(3) - S 1)2 (x)] (133) 

= x{Si )2 (l - x) - Si )2 (x) 



X 



ln 2 (x) 
1 — x 



X 



'ln 3 (l-x)' 
1 -x 



+ 



Q ln(x)ln 2 (l-x) 
xj + 1 — X 



ln 2 (x) ln(l -x) 
xj + 1 — X 



+ 



X - ln(x) ln(l - x) [ln(l - x) - ln(x)] + C(3)| (134) 



x 



+ ln(x) ln(l - x) + ^ ln 2 (x) + (1 - x) ln(l - x) 
+x ln(x) + xLi 2 (l — x) 



= x 



2S li2 (l-x)--ln 3 (x) 



+ ln 2 (x) 



+2[ln(x) + 1 - x] 
51n 4 (l-x) , . > 3 (l-x) 
A-T^x HX) 1 



+ 6 C (3)M1^)_ 6C (4) 

X Ob 



1 - x 
1 



1 -x 



1 - X 



-2S 2 , 2 (x) - 2S 2 , 2 (1 - x) + 



3C(2) 



C(4) 



ln 2 (l -x) 
1 -x 



(135) 



(136) 



(137) 



-Li 2 (l - x) + 6 ln(l - x)S 1>2 (l - x) + 4 ln(x)Si, 2 (x) 
-2 ln 2 (x) ln 2 (l - x) + ln(x) ln 3 (l - x) 

-2C(2)ln(x) ln(l-x) 

41n(l -x)Si, 2 (l -x) -Li 2 (l -x) 2 
+2 ln(x)[S 1;2 (l - x) + S li2 (x) - C(3)] (138) 
- ln 2 (x)Li 2 (x) - \ ln 3 (x) ln(l - x) 
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1 \ ln 3 (x) 



1 — X J _|_ 1 — X 



1 — X/ + 1 — X 



1 



8 Li 3 (l-,) 

.1 — X J + 1 — X 

1 \ _ ln(x)Li 2 (l — x) 



x 



1 — X 



J_\ s Li 3 (x)-C(3) 
1 — xj + 1 — X 



1 — X 



1-x 



+ 



1 



1 — X/ + 

1 



1 — X 



+ 



g)Li 3 (l -x) 
xLi 3 (l — x) 



1-x 



-S li3 (l-x) + C(4)-Li 4 (x) 
1 _ . o, , 1 



+C(3)ln(x) + -C(2)ln 2 (x)--ln 4 (x) 



2S 2>2 (l-x)-2S li3 (l-x) 



(139) 



1-x 
1 



2 Li^(l-x)+ln(^)S l2 (l-,) 



In ( ^— ^ ) Li 3 (l - x) - S 2)2 (l - x) 



1 — x 
1 

1 — X 



X 



(140) 
(141) 



2S 2 , 2 (x)-2S 2 , 2 (l-x)+4S 1;3 (l-x) 



+_ Li 2 (1 _ x) + ln(x)[S li2 (l - x) - S 1)2 (x) - C(3)] 
~ ln 2 (x)Li 2 (l - x) - 2 ln(l - x)Si, 2 (l - x) 



C(4) 



-2S 2 , 2 (x) - 2Si, 3 (l - x) 



1 — x 

-ln(x)[S li2 (l-x)-2C(3)] 

- ln(l - x)[S 1)2 (l - x) - C(2) ln(x)] + ±C(4) 

+iLi 2 (l-x) 2 + iln 3 (x) ln(l-x) 



g,S 1)2 (x) = S 2 , 2 (x)-3S 1 , 3 (x)-ln(x)S 1 , 2 (x) 

+C(2)Li 2 (x) + C(3)ln(l-x) + ic(4) 
xSi )2 (x) = x|s 2j2 (x) - 3Si )3 (x) - ln(x)Si )2 (x) 

C(2)Li 2 (x)+C(3)ln(l-x) + ic(4)| 



(142) 
(143) 



(144) 



+C(3)-zS 1;2 (x 
-(1-x) 



l-ln(l-x) + -ln 2 (l-x) 



-S 2 , 2 (l - x) - Li 4 (l - x) + In 

£ j-S 2j2 (l - x) - Li 4 (l - x) 
'1 — X s 



1 — X 



X 



(145) 

Li 3 (l - x) (146) 



+ ln 



x 



Li 3 (l-x) 



+xln(x) + (1 - x)[l - Li 2 (l - x) + Li 3 (l - x)\ (147) 
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ln(l — x)Li 2 (l — x) 



2S 2 , 2 (x) - S 1)3 (x) + 



C(4) 



(148) 



1 



-3 ln(l - x)Si >2 (l -x)-2 ln(x)S li2 (x) + |Li 2 (l - x) 2 



+ 



■ ln 2 (l - x) - ln(x) ln(l - x) - C(2) 



Li 2 (l-x) 



1 1 
-- ln(x) ln 3 (l - x) + - ln 2 (x) ln 2 (l - x) 

vJ ^ 



xln(l -x)Li 2 (l -x) = x^2S 2 , 2 (x) - Si, 3 (x) + 



C(4) 



(149) 



ln(x)Li 2 (l — x) 



xln(x)Li 2 (l — x) 



-3 ln(l - x)Si, 2 (l -x)-2 ln(x)Si, 2 (x) + -Li 2 (l - xf 
+ iln 2 (l-x)-ln(x)ln(l-x)-C(2) Li 2 (l - x) 

-\ ln(x) ln 3 (l -x) + l ln 2 (x) ln 2 (l - x)\ 
6 2 j 

+ (1 - x) [ln(l - x) - 1] Li 2 (l - x) - Li 2 (l - x) 

+3(1 — x) + 2x ln(x) — [1 — x + x ln(x)] ln(l — x) 

2S 2 , 2 (x) - 2S 2>2 (1 - x) + 4Si )3 (l -x) + Li 2 (l - x) 

+ ln(x)[S 1)2 (l-x)-S li2 (x)-C(3)] 

~ln 2 (x)Li 2 (l-x) 

-21n(l-x)S 1)2 (l-x)-ic(4) (150) 



x- 



2S 2 , 2 (x)-2S 2 , 2 (l-x)+4S li3 (l-x) 



+ ln(x)[S 1)2 (l-x)-S 1;2 (x)-C(3)] 
+Li 2 (l - x) - - ln 2 (x)Li 2 (l - x) 



1-x 



Si, 2 (l-x) 



-21n(l-x)S 1)2 (l-x)--C(4)| 

+ [1 — x + xln(x)](3 — Li 2 (l — x)) — xln 2 (x) 
-2S 1)2 (l-x) 

-2S 1)3 (1 -x) + S 2>2 (1 - x) - ^Li 2 (l - x) 
'1 — x 



+ ln (— *) S li2 (l-x) 



®xS 1;2 (l-x) = 

1 — X/ + 



X 



-2S li3 (l-x) + S 2 , 2 (l-x)--Li 2 (l-x) 



+ ln('^) S li2 (l-x) 



+(l-x)[Si, 2 (l -x) - 1] -xln(x) 



1 - -ln(x) 



(151) 



(152) 



(153) 



22 



,1 -x 



+ 



(- ) ®xLi 3 (x) = x 

VI — x/ + 



Li 3 (x) = -2S 2 , 2 (x)-2S li3 (l-x)-ln(x)S li2 (l-x) 

-C(2)Li 2 (a;) + i ln 3 (x) ln(l - x) + 3C(4) 
o 

+C(3)[ln(x)+ln(l-x)] 

2S 2 , 2 (x) - 2Si, 3 (l -x)- ln(x)S li2 (l - x) 

<(2)Li 2 (x) + - ln 3 (x) ln(l - x) + 3C(4) 



— a; 



+ 



- x 



x ln(x)Li 2 (x) 



— ^— ) <g>ln 3 (l-x) 
1 — x/ + 



— x 



■ xln 3 (l — x) 



<S> ln(x)Li 2 (x) = 



+C(3)[ln(x)+ln(l-x)] 

+C(3)-C(2) + z[Li 2 (z)-Li 3 (z)] 
+(1 -x)[l -ln(l - a;)] 

-2S 2 , 2 (x)-4S 1)3 (l-x) 

+ ln(x)[3C(3)-S 1 , 2 (x)-3S 1 , 2 (l-x)] 

+iln 2 (x)Li 2 (x) + 3C(4)-C(2)Li 2 (x) 

+^ln 3 (x) ln(l -x) 

xj-2S 2 , 2 (x)-4S li3 (l-x) 

+ ln(x)[3C(3)-S 1 , 2 (x)-3S 1 , 2 (l-x)] 

+iln 2 (x)Li 2 (x) + 3C(4)-C(2)Li 2 (x) 

+^ln 3 (x) ln(l - x)j + [1 + a; — xln(x)]Li 2 (x) 

+ (1 - x)(2 - ln(x))[l - ln(l - x)] 
-2C(2) + 1 -x + ln(x) 

ln 4 (l - x) - ln(x) ln 3 (l - x) - 3({2) ln 2 (l - x) 
+6C(3) ln(l - x) - 6C(4) 

x ln 4 (l - x) - ln(x) ln 3 (l - x) - 3((2) ln 2 (l - ; 



+6C(3) ln(l - x) - 6C(4) 



-6(1 -x) 
1 



1 



l-ln(l-x) + -ln 2 (l-x) 
2 



® ln(x) ln 2 (l — x) = 



-ln 3 (l -x) 

6 v 1 



-2S 2 , 2 (x) - 2S 2>2 (1 - x) - 2S 1)3 (x) + -C(2) 
-Li 2 (l - x) + 6 ln(l - x)S li2 (l - x) 



23 



-41n(x)[Li 3 (l-x)-C(3)] 

+2 ln(x) ln(l - x) [2Li 2 (l - x) - C(2)] 
2 

+ -ln(x) ln 3 (l - x) 
3 

xln(x)ln 2 (l-x) = x{-2S 2 , 2 (x)-2S 2 , 2 (l-x)-2S 1)3 (x) + ^C(2) 

-Li 2 (l - x) + 6 ln(l - x)Si )2 (l - x) 

-41n(x)[Li 3 (l-x)-C(3)] 

+2 ln(x) ln(l - x) [2Li 2 (l - x) - C(2)] 

+ ln(x) ln 3 (l - x) | - 2x ln(x) 

+ (1 - x) ln(x) ln(l - x)[ln(l - x) - 2] 
-(1 - x)[ln(l - xf - 41n(l - x) + 6] 
+2[C(2) - Li 2 (x)] + 2[C(3) - S 1>2 (x)] 

= 2S 2)2 (1 -x) + 2[ln(x) + ln(l - x)]S 1)2 (l - x) 

+21n(x)[S li2 (x)-C(3)]-Li 2 (l-x) 2 

+ ln 2 (x) [Li 2 (l - x) - C(2)] - \ ln 3 (x) ln(l - x) 



^ <g>ln 2 (x)ln(l -x) 



xln 2 (x)ln(l -x) = x|2S 2i2 (l -x) +2[ln(x) +ln(l -x)]Si )2 (l -x) 

+21n(x)[S li2 (x)-C(3)]-Li 2 (l-x) 2 
+ ln 2 (x) [Li 2 (l - x) - C(2)] - i ln 3 (x) ln(l - x) j 
+2S 1;2 (l-x) + 2Li 2 (l-x) 

-x[ln 2 (x) - 2 ln(x)] ln(l - x) + 2(1 - x) ln(l - a 
-4xln(x) +xln 2 (x) + 6(x - 1) 



1 -x 



ln 3 (x) 



ln(x) ln 2 (l — x) 
1 — x 



6 



C(4) - Li 4 (x) + C(3) ln(x) + -C(2) ln 2 (x) 



-^ln 4 (x) 



f — 3—) <g>xln 3 (x) = 6x 
VI — x/ + 



C(4)-Li 4 (x) + C(3)ln(x) + -C(2)ln 2 (x) 



+ 1 - ln(x) + ]- ln 2 (x) - i ln 3 (x) 

-6--ln 4 (x) 
4 

6[S li3 (x)-C(4)] + ^ln 4 (l-x) 



= 2 



S 2 , 2 (x) + S li3 (x) - ln(x)S 1;2 (x) - -C(4) 
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+-ln(x) ln 3 (l - x) 



ln 2 (x) ln(l - x) 



1-x 



, In 3 (:r) 
1 <g> — — 

1 — X 

1 — X 

1 — X 
1 ln(x)Li 2 (l -ap 

1 — X 

10 LijM-C(3) 

1 — X 



l®Si, 2 (x) = 

1<8xSi )2 (x) = 

1 <8 Li 3 (l - x) = 

1 <8 xLi 3 (l — x) = 

1 ® ln(l - x)Li 2 (l - x) = 

1 <S> xln(l — x)Li 2 (l — x) = 

1 (8) ln(x)Li 2 (l - x) = 

1 <8 xln(x)Li 2 (l — x) = 



-2 



S 2 , 2 (l-x) + S 1;3 (l-x) 
1 



(167) 
(168) 



ln(l - x)Si, 2 (l - x) + | ln 3 (x) ln(l - x) 
-lln 4 (x)-6Si )3 (l-a;) 



S 2)2 (l - x) - ln(x)Si, 2 (l -x)- 3S 1)3 (1 - x) 

— -Li^Cl - x) - ln(x)Li 3 (l - x) + Li 4 (l - x) 

-|Lil(l " a) " ^ ln 2 (x)Li 2 (l - x) + 3S lj3 (l - 

C(4) - Li 4 (x) + C(3) ln(x) + ln(l - x)[Li 3 (x) - 

-i[C(2) 2 -Li 2 (x)] 

|C(4) - S 2 , 2 (x) 
C(3)-xS li2 (x) 

-(1-x) l-ln(l-x) + iln 2 (l-x) 

— ln(x)Li 3 (l — x) — 2^^2(1 ~~ x ) 

xln(x) + (1 - x)[l - Li 2 (l - x) + Li 3 (l - x)] 
2S 2j2 (l - x) - iLi 2 (l - x) - 2 ln(l - x)S 1;2 (l 

— ln(x) ln(l — x)Li 2 (l — x) 

-Li 2 (l - x) + (1 - x)[ln(l - x) - l]Li 2 (l - x) 
+ (1 - x) [3 - ln(l -x)]+x ln(x) [2 - ln(l - x) 

-iln 2 (x)Li 2 (l-x) + 3S li3 (l-x) 

3(1 — x) + x ln(x) [3 — ln(x)] 

-x[ln(x) - l]Li 2 (l - x) - 2Si, 2 (l -x)- Li 2 (l 



(169) 
(170) 
(171) 
x) (172) 

C(3)] 
(173) 
(174) 

(175) 

(176) 
(177) 



— x) 



(178) 

(179) 
(180) 



— x) 



(181) 



l®Si )2 (l-x) 


= -ln(x)S 1)2 (l-x)- 


3Si, 3 (l 


-x) 


(182) 


1 ®xSi i2 (l -x) 


= (l-x)[S 1)2 (l-x)- 


1) -x 


ln(x) ln 2 (x) 


(183) 


1 <g> Li 3 (x) 


= C(4)-Li 4 (x) 






(184) 


1 ® xLi 3 (x) 


= (1 -x)[l -ln(l -x 


)] +x[Li 2 (x) -Li 3 (x)] 






+C(3)-C(2) 






(185) 


1 <8 ln(x)Li 2 (x) 


= Li 4 (x) — ln(x)Li 3 (x) 


- C(4) 




(186) 


1 <8 xln(x)Li 2 (x) 


= [3-21n(l-x)](l - 


x) - Li 2 (l - x) - C(2) 
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1 <g>ln 3 (l -x) 
1 <g>xln 3 (l - x) 

1 ® ln(x) ln 2 (l - x) 
1 ® xln(x) ln 2 (l — x) 



1 <g> ln 2 (x) ln(l -x) 
1 ® xln 2 (x) ln(l — x) 



+x {Li 2 (x) [1 - ln(x)] + ln(x) [1 - ln(l - x)]} 
= -6[C(4) - S 1;3 (x)] 
= -(l-x)[6-61n(l-x)+31n 2 (l-x) 

-ln 3 (l -x)] 

= -2^C(4)-S 2 , 2 (x)+ln(x)S li2 (x) 

= 2[C(3) - S li2 (x)] + 2[C(2) - Li 2 (x)] - 2 ln(x) 
-(1 - x)[6 - 41n(l - x) + ln 2 (l - x)] 
+ (1 -x) ln(x)[2-21n(l - x) + ln 2 (l - x)] 

= -2[S li3 (l-x) + iln 3 (x)ln(l-x)] 



x 



X 



'ln 3 (l-x) 
1 — x 



ln(x) ln 2 (l — x) 
1 -x 



x ® 



ln 2 (x) ln(l — x) 
1 — x 



x ® 



ln 3 (x) 
1 -x 



x ® 



Si, 2 (l-x) 
1 -x 



S 1;2 (l-x)+Li 2 (l-x) 



+xln(x)< -2 + 



2 - ln(x) 



ln(l-x) - 1 



+2(1 -x) [ln(l-x)-3] 

6x [ Sli3 (x) + S li2 (x) - C(4) - C(3)] 



+ |ln 4 (l-x) + (l-x) ln 3 (l-x) 

2x S 2 , 2 (x) + S 1;3 (x) - ln(x)S li2 (x) - ^C(4) 

-S 1)2 (l - x) + S li2 (x) + 1 ln(x) ln 3 (l - x) 

~ ln(x) ln(l - x) [ln(l - x) - ln(x)] - C(3) 

+ (1 - x) ln 2 (l - x) - 2x [Li 2 (x) - C(2)] 
+ ln(x) ln 2 (l - x) 



= -2x 



(187) 
(188) 

(189) 
(190) 



(191) 
(192) 



(193) 



(194) 



(195) 



S 2)2 (l -x) + S li3 (l -x)- ln(l - x)S li2 (l 

1 



+Si, 2 (l - x) - Li 2 (l - x) + | ln 3 (x) ln(l - x 

— - ln 3 (x) — — ln 2 (x) + ln(l — x) — ln(x) 

+2[1 + ln(x)] ln(l - x) + ln 2 (x) ln(l - x) 
"I 



x) 

(196) 



— x 



+6 



ln 4 (x) + 6S li3 (l-x) 



1 - x + ln(x) + ^ ln 2 (x) + i ln 3 (x) 



X 



S2,2(l - x) - ln(x)Si, 2 (l - x) - 3Si ;3 (l - x) 



x 



+ (l-x)S li2 (l-x) + ^ln 3 (x) 



(197) 



(198) 
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X 



Li 3 (l-x) 
1-x 



X ® 



ln(x)Li 2 (l — x) 



1 — x 



x ® 



Li 3 (x)~C(3) 
1-x 



x ® Si )2 (x) 

x ® xSi )2 (x) 
a; (g) Li 3 (l — x) 

x ® xLi 3 (l — x) 
x <g> ln(l — x)Li 2 (l — x) 



x <g> ln(x)Li 2 (l — x) 

x <g> xln(x)Li 2 (l — x) 

x ® Si )2 (l - x) 
x ® xSi j2 (l — x) 



— x 



-Li 2 (l - x) + ln(x)Li 3 (l - x) - Li 4 (l - x) 



+ (1 — x)Li 3 (l — x) + x ln(x)Li 2 (l — x) 
+2xSi, 2 (l -x) 



x 



_I Li 2 (1 -x)- 1 - ln 2 (x)Li 2 (l - x) + 3S 1)3 (1 



X, 3 s X 



+ [1 + ln(x)]Li 2 (l - x) - ^ ln 3 (x) - ^ ln 2 (x) 
+2xSi, 2 (l - x) - xLi 2 (l - x) 

= Li 3 (x) + Li 2 (x) - C(3) - xC(2) - x ln(x) 
— (1 — x) ln(l — x) 

+x|c(4)-Li 4 (x)+C(3)ln(x) 
+ ln(l-x)[Li 3 (x)-C(3)] 
-i[C(2) 2 -Li 2 (x)]j 

= S 1)2 (x) - xC(3) + \{l-x) ln 2 (l - x) + x[C(2) 

-Li 2 (x)] 

= x gc(4) - S 2j2 (x) 

= 2xS li2 (l - x) + (1 - x)Li 3 (l - x) 
+x ln(x)Li 2 (l — x) 

= —x ln(x)Li 3 (l — x) + -Li 2 (l — x) 

= 3xS li2 (l-x)-^ln 2 (x)ln(l-x) 

+ [x ln(x) + (1 - x) ln(l - x)] Li 2 (l - x) 



x <g> x ln(l — x)Li 2 (l — x) = x 



2S 2 , 2 (l-x)--Li 2 (l-x) 



-21n(l -x)Si )2 (l -x) 
-ln(x) ln(l -x)Li 2 (l - x) 

2xSi, 2 (l - x) + [1 - x + ln(x)]Li 2 (l - x) 

--ln 3 (x) - -ln 2 (x) 
3 2 



— X 



^ln 2 (x)Li 2 (l-x)-3S li3 (l-x) 



(l-x)S li2 (l-x) + -ln 3 (x) 
o 



= x [— ln(x)Si j2 (l — x) — 3Si j3 (l — x)] 
x®Li 3 (x) = Li 3 (x) + Li 2 (x) -x[C(3) + C(2)] -xln(x) 



(199) 

x) 

(200) 



(201) 

(202) 
(203) 

(204) 
(205) 

(206) 



(207) 

(208) 
(209) 

(210) 
(211) 
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x <S> xLi 3 (x) 
x <S> ln(x)Li 2 (x) 



x <S> x ln(x)Li 2 (x) 
x <g) ln 3 (l - x) 
x <g> xln 3 (l — x) 



x ® ln(x) ln 2 (l — 



x) 



— (1 — x) ln(l — x) 
x[C(4)-Li 4 (x)] 

-xLi 2 (l - x) + [1 + ln(x)]Li 2 (x) - <(2) 
-ln(l -x)[ln(x) + 2] 

— 2a; ln(x) — ln(l — x) + — ln 2 (x) 

x [Li 4 (x) - ln(x)Li 3 (x) - C(4)] 
6x[S li2 (x)-C(3)] + (l-x) ln 3 (l- 
-6x [C(4) - S 1)3 (x)] 

[1 - x + ln(x)] ln 2 (l - x) + 2x[C(2) - Li 2 
-2x C(3) - Li 3 (x) + ln(x)Li 2 (x) + Li 3 (l - x) 



ln(l - x)Li 2 (l - x) 



■ x) 



(212) 
(213) 



(214) 

(215) 
(216) 
(217) 



x ® xln(x) ln 2 (l — x) 
x <g> ln 2 (x) ln(l - x) 



-2x 



-C(4)-S 2l2 (x)+ln(x)S 1|2 (x) 



(218) 
(219) 



2x [Li 2 (l - x) - S 1)2 (l - x)} 

2[1 - x + ln(x)] + ln 2 (x)] ln(l - x) 
1 



+ 

+x ln(x) 



2 + ln(x) + ^ln 2 (x) 



x ® x ln 2 (x) ln(l — x) = 



-2x[S li3 (l-x) + -ln 3 (x) ln(l-x)] 



(220) 
(221) 



For some applications one wishes to give explicit account for soft gluon or soft photon re- 
summation using x-space formulae. Here one has to perform convolutions of the type 



"ln fc (l -x)~ 




V(l -x)~ 


1 -x 


+ 


1 — X 



D k (x)<g>Di(x), k,l>0 



(222) 



The Mellin convolution of two (...) + -distributions is again a (...) + -distribution. The convolutions 
Dk <8> Dq, k e [0, 3] are given in Eqs. (69,72,104) and (131) above and can easily be extended to 
general values of k and / using relations given in Ref. [21]. The Mellin transform of D k (x) reads 



M[D k (x)](N) 



-lnk-iy.s^N-i). 

k 



(223) 



The harmonic sum in (223) obeys a determinant-representation, Eq. (158) [21], and is a polyno- 
mial of only single harmonic sums. One obtains the Mellin transform of [ln(x) ln^(l— x) / (1— x)]+ 
by the relation 



M 



ln(x 



ln fe (l -x) 



x 



(N) 



_d_ 

dN 



-l) k (k-l)\S 1 ,..., 1 (N-l). 



(224) 



Since the single harmonic sums can be represented by Euler's ^-function and its derivatives the 
differentiation for N is performed easily. Furthermore, one may re-express the result again in a 
polynomial of single harmonic sums. 
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9 Appendix B: Mellin Transforms 

Some of the convolutions calculated in the previous section were determined using harmonic 
sums and Mellin transformations of specific functions, such as Nielsen-integrals. The Mellin 
transform of a convolution is given by 



M{[A® B](x)}(N) = M{[A](x)}(N) ■ M {[B](x)} (N) . 
Below we list Mellin transforms which were used and were not contained in Ref. [21]. 



(225) 



M[ln p (l -x)](N) 
M[ln 4 (l - x)](N) 



M[S ltP (x)](N) 

M[S hp (l - x)](N) 
M[S 1)3 (x)](7V) 

M[Si )3 (l-a;)](iV) 

M[S 1A {x)](N) 



M 



M[S 1A (l-x)](N) 

M[S 2 , 2 (x)](N) 

M[S 2t2 (l-x)](N) 
S 2 , 2 (l - x 



1 — x 



M 



Ufa) (N) 



±(-iy p \ Sl _i(N) (226) 
p 

^ [st(N) + 6Sf(N)S 2 (N) + 3S 2 2 (N) + 8S1(N)S 3 (N) 
+6S 4 (N)\ 



1 

N 



C(f)-/i 



-„ [S P (N) - Cip)} 



Sf(N) + 3S 1 (N)S 2 (N) + 2S 3 (N) 



N 



[S 4 (N) - C(4)] 



(227) 
(228) 

(229) 
(230) 
(231) 



11 1 

iV ~ 24/V 



Sl(N) + 6Sf(N)S 2 (N) + 8S 1 (N)S 3 (N) 



3S%(N) + 6S 4 (N) 



C(4) C(3) , S 2 (N) + S 2 (N) 



+ 



2N 3 



AN N 2 

7C( 4 ) + Si,s(N) - CmSi(N) 



1 

N 



(N) = 2C(5) - C(2)C(3) - -^WSiiN - 1) - S hl , 3 (N - 1) 



+C(3)5i,i(JV 
AM 2 SV ; 



N 



N 



N 



+ 



N 3 



S*(N) + S 2 (N) 



(232) 
(233) 
(234) 
(235) 

(236) 

(237) 



M [ln(l — x)Li 3 (x)] (N) = 



1 

TV 



1 

2 

S 2 ,i(N) 



C(2) 2 + ^P + C(2) 



S 1 (N) S 2 {N) + S 2 {N) 



N 



N ' N N 2 

C(3)S 1 (N) + C(2)S 2 (N) - S^(N) 



(238) 
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M[\n(x)U 3 (x)}(N) = -± 



' 3C(2) | S 2 (N) | 3S 1 (N) 



N 



N 



N 2 



M[\n(x)S lt2 (l-x)} (N) 



M[ln(x)S 1;2 (x)] (AO 



^{3^4(iV)-3C(4) + l[53(iV)-C(3)]} 
-^U(N)S 2 (N) + S 3 (N) - C(2)5x(JV) - 2C(3) 



(239) 
(240) 



+ 



1 r 

N 



S 2 (N) + S 2 (N) 



M [Li 4 (x)} (N) 



C(4) C(3) + C(2) S^N) 



N 



N 2 N 3 



N 4 



(241) 
(242) 



M 



M 



M 



\n 2 (x)Uo(l -x) 



S 2 , 2 (x)-C(4)/4 
1 -x 

S 1)3 (x)-C(4) : 



(TV) 



i|6[C(4)-5 4 (JV)] + ^[C(3)-5 s (JV)] 



M 



1 — x 



(TV) 

(AO 
(AO 



+_[C(2)-S 2 (JV)]j 
-5 3il)1 (JV-l) + C(3)5 2 (JV-l) 



5 2>1>1>1 (JV - 1) - 4C(5) 



C(5) 



5 M (iV - 1) - C(4)5i(JV - 1) + 2C(5) - C(2)C(3) 



(243) 
(244) 
(245) 
(246) 



10 Appendix C: Representation of Sums 



To obtain the inverse Mellin transforms of more complicated convolutions back in x-space prod- 
ucts of harmonic sums and powers of 1/N have to be expressed as Mellin transforms. A series 
of representations has been given in [21] before. Here we list additional representations which 
were used in the present calculation. 



giCjv) 

N 
Si(N) 

N 2 
Si(N) 

N 3 

N 4 

S 2 (N) 

N 
S 2 (N) 

N 2 
S 2 (N) 

N 3 



-M[ln(l -x)](N) 

-M[Li 2 (z)-C(2)](A0 

M[Li 3 (x)-C(2)ln(x)-C(3)](7V) 



— M 



LU(x)-^-\n 2 (x)-((3)Hx)-((4) 



(N) 



-M[U 2 (l-x)-((2)](N) 

M[2S 1)2 (1 -x) + ln(x)[Li 2 (l -x)- C(2)]](iV) 

1 



= — M 



3S 1)3 (1 -x) + 2 ln(x)S li2 (l -x) + - ln 2 (a;)[Li 2 (l - x) - C(2)] 



(247) 
(248) 
(249) 
(250) 

(251) 
(252) 

(AO 
(253) 
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SsjN) 
N 

N 2 
S*(N) 
N 

sm 

N 

sm 

N 2 

sm 

N 3 



sm 

N 

S 2 (N)S 1 (N) 
N 

S 2 (N)S 1 (N) 
N 2 

S 3 (N)S 1 (N) 
N 

s 2 (N)sm 

N 



-M[S 1)2 (l-x)-C(3)](A0 

M[3Si, 3 (l -x) + In(x)Si >2 (l - x) - C(3) ln(x)](iV) 

- M[8 lfl (l - x) - C(4)](N) 

M[ln 2 (l - x) + Li 2 (l - x) - ((2)}(N) 



(254) 
(255) 
(256) 
(257) 



M[2S li2 (x) + 2Si, 2 (l - x) + ln(x)[Li 2 (l - x) - C(2)] - 2C(3)](iV) (258) 
C(4) 



M 

+ M 



2S 2 , 2 (x) - 



-2C(3)ln(x) 



(N) 



3S 1>3 (1 -x) + 2 ln(x)S li2 (l ~x) + - ln 2 (x)[Li 2 (l - x) - C(2)] 



M [4S 22 (x) + S 13 (l - x) - 4 ln(x)S 1;2 (x) - 2 ln(x) ln(l - x)Li 2 (x)] 
+ M [4 ln(l - x)Li 3 (x) + Li 2 (x) - 4£(3) ln(l - x) - C(4)" 

M [S 1)2 (l - x) - Li 3 (l - x) + ln(l - x)[Li 2 (l - x) - C(2)]] (iV) 

M [ln(x)[S li2 (x) - S 1)2 (l - x)\ - 3S 1>3 (1 - x)] (iV) 
- M [2S 2)2 (x) + C(2)Li 2 (x) - C(3) ln(x) - 3((4)] (N) 



= M 



S 2 , 2 (l-x) + S li3 (l-x)--Li 2 (l 



x 



(N) 



+ M[[ln(l-x)[S 1)2 (l-x)-C(3)]](7V) 

M [-4S 22 (x) - 2S 22 (1 - x) + 2Si 3 (x) - S 13 (l - x) + 41n(x)S li2 (x)] 



(iV) 
(259) 

(260) 
(261) 

(262) 
(263) 



-2 ln(l - x)S li2 (l - x) - 4 ln(l - x)Li 3 (x) + Li 2 (l - x) - Li 2 (x 



+ M 
+ M 

+ M[4C(3)ln(l-x)-C(4)] 



2 ln(x) ln(l - x)Li 2 (x) + \ ln(x) ln 3 (l - x) + C(2) ln 2 (l - x) 



(264) 
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Figure 1: The non-singlet radiation function D^ s (x,Q 2 ) to 0((aL) 5 ) as a function of x for Q = 
lOGeV, lOOGeV and ITeV. 
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Figure 2: The effect of the resummed contributions beyond 0((aL) 5 ) compared to the contributions 

up to 0((aL) 5 ) for Dns(x, Q 2 ) as a function of x and Q. 
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Figure 3: Relative contribution of the first order non-singlet radiator compared to all terms up to 
0((aL) 5 ) as a function of x and Q and relative correction of all contributions up to 0((aL) 4 ) in 
comparison to the terms up to 0((aL) 5 ). 




Figure 4: The polarized singlet contributions Dij as a function of a: and Q in %. Here D\\ denotes 
the pure singlet term, which has to be added to the non-singlet contribution. 
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Figure 5: Relative contribution of the first order singlet radiators in the terms to 0((aLf 
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Figure 6: Relative contribution of all terms of the singlet radiators Dij to the 4th order in aL if 
compared to all terms to 5th order as a function of x and Q. 
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Figure 7: The polarized singlet contributions Dij(x,Q 2 ) corresponding to the resummation of the 
0(aln 2 (x)) terms as a function of x and Q in % starting with 0{a 2 ). Here D\\ contains also the 
non-singlet contribution. 




Figure 8: The QED radiative initial state correction to the differential cross section of polarized 
deep-inelastic lepton-proton scattering including the 0((aL) 5 ) corrections as a function of x and y 
in a fixed target experiment (E e = 27.5 GeV). Left : the correction factor. Right : the contributions 
due to non-leading orders. 
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Figure 9: The QED radiative initial state correction to the differential cross section of polarized 
deep-inelastic lepton-proton scattering including the 0((aL) 5 ) corrections as a function of a; and y 
for a collider experiment at S — ITeV 2 . Left : the correction factor. Right : the contributions due 
to non-leading orders. 
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